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^ ; Abstract 

Quantum Defect theory is a well established theoretical concept in modern spectroscopy. We 



show that this approach is useful in electron impact ionization problems where state of the art 
theoretical methods are presently restricted mostly to simple atomic targets. For the well docu- 
mented Argon ionization case in equal energy sharing geometry the approach suggested leads to 



X 

■ significant improvements compared to previous calculations. 
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I. INTRODUCTION 



Quantum Defect (QD) information is widely exploited in modern spectroscopy, to 
characterize Rydberg states and in the calculation of the photoionization cross sections of 
various atomic and molecular species [|^, 0. 

In the present work, we show that QD information, used within the Distorted Wave 
Born Approximation (DWBA) framework, might also prove useful for the description 
of ionization processes by presenting an alternative way to account for the short range 
interactions (static and exchange) in the calculation of the final state continuum distorted 
waves. The range of validity of this approach reaches beyond that |^ of the commonly 
used Furness-McCarthy local exchange approximation [Q. Compared to the determination 
of the Hartree-Fock non-local operator [H, which becomes rapidly a tedious task as the size 
of the target increases, our method allows for a target-independent procedure which can be 
readily applied to much larger atomic or molecular systems. 

The parameters of the Green-Sellin-Zachor parametric form of the electron-ion 
potential are optimized in order to reproduce the QD using the canonical function method 
1^. The canonical function method (CFM) is a powerful means for solving accurately the 
Radial Schrodinger Equation (RSE). The quantum eigenvalues of the RSE are obtained to 
any desired accuracy without having to evaluate the eigenf unctions. Several parametric 
potentials are discussed extensively in the literature and we already used several functional 
forms adapted to different atomic systems. We studied recently |^ the QD of some 
rare gases with the Klapisch parametric potential and found that in some cases it 
was very difficult to optimize parameters that provide an accurate representation of the 
experimental QD. We believe, the Green-Sellin-Zachor @ is better suited to our present 
study as the parameter space is small (two-dimensional) which allows for an efficient search 
of the optimized parameters. 

The electron-ion potentials obtained for each Rydberg series are further modified to 
account classically for the electron-electron interaction in the final state. Calculations per- 
formed within the DWBA framework for the ionization of argon in the equal energy sharing 
geometry = — k^, which is reasonably well documented both theoretically and experi- 
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mentally, clearly validate our approach and shows significant improvements over previous 
treatments. 

This paper is organized as follows: Section 2 is devoted to the presentation of the CFM, 
section 3 describes the optimization method we use based on the CFM. In section 4 we 
describe the results whereas section 5 contains our conclusion and discussion of the results. 



II. THE CANONICAL FUNCTION METHOD 

The canonical function method (CFM) |]^ is a powerful means for solving the Radial 
Schrodinger Equation (RSE). The mathematical difficulty of the RSE lies in the fact it is 
a singular boundary value problem. The CFM turns it into a regular initial value problem 
and allows the full determination of the spectrum of the Schrodinger operator bypassing 
the evaluation of the eigenf unctions. 



The Canonical Function Method (CFM) developed by Kobeissi |]T0[ and his coworkers 
to integrate the RSE, consists in writing the general solution as a function of the radial 
distance r in terms of two basis functions a{E] r) and (3{E] r) for some energy E. 



Generally, one avoids using the wavefunction but if one needs it, the functions a{E] r) 
and l3{E; r) are used to determine the wavefunction at any energy with the expression: 

^{E; r) = ^{E; ro)a{E; r) + ^'{E; r^)(3{E- r) (1) 

where '^{E\ tq) and '^'{E] tq) are the wavefunction and its derivative at some initial distance 



At the selected distance rg a well defined set of initial conditions are satisfied by 
the canonical functions and their derivatives ie: a(i?;ro) = 1 with a'(i?;ro) = and 
/?(£;; ro) = with /3'(E;ro) = l. 

The method of solving the RSE is to proceed from tq simultaneously towards the origin 
(r = 0) and towards oo. During the integration, the ratio of the canonical functions is 
monitored until saturation signaling the stability of the eigenvalue spectrum. 
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The saturation of the ratio with r yields a position independent eigenvalue function 



F{E). The latter is defined with the help of two associated energy functions: 

a{E-r) 



and: 



1+{E) = linir^+c 



1~{E) = lirrir- 



(3{E;r) 

a{E;r) 
WW) 



(2) 



(3) 



The eigenvalue function is defined in terms of: 



F{E) = l^{E)~l_{E) (4) 

Its zeroes yield the spectrum of the RSE. An example of behaviour of F{E) is displayed 
in Figure 1. The eigenfunctions may be obtained for any E = Ek where E^ is a zero of F{E). 
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FIG. 1: Typical behavior of the eigenvalue function with energy. The vertical lines indicate the 
eigenvalue position. 

The speed and accuracy of the method have been tested and compared to standard 
integration algorithms in a variety of cases and for a wide of range of potentials. In addition 
the method has been generalised to the phase shift estimation as well to the continuum 
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coupled channel problem 
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In potential optimisation problems the method can only be of great value because of 
its speed, accuracy and ease of programming. Accordingly, one has to pick an appropriate 
optimisation method and use the CFM as a subroutine for a given parametric potential. 



III. OPTIMISATION PROCEDURE 



An extensive review of the applications of model potentials has been given by Hibbert 



11| and Aymar et al. H. The functional form suggested by Green et al. p| is given by: 



\/(r) = -(2/r)[(Z-lV(r) + l], 
with a;(r) = l/[ei(exp(r/e2) - 1) + 1] (5) 

where ei and €2 are parameters that are determined by the optimisation procedure. 



The optimisation problem, at hand, is over-determined since the experimental set of 
energy levels might consist of tens of values whereas the potential depends only on two 
numbers, namely ei and 62. This over-determination allows us to use several criteria for 
the optimisation procedure and later on select the best one that achieves results closest to 
experiment. 



The optimization procedure consists of defining an objective function and finding its 
minimum in the two-dimensional parameter space ei, €2- The objective function is a 
quadratic consisting of the difference between some picked levels and those produced by the 
parametric potential through the solution of the RSE with the CFM. 



We adopted several strategies based on the following observations: The quantum defect 
value is not stable for the low levels but tends to reach a stable value when the energy 
increases. Despite the stability of the quantum defect for the higher levels, the experimental 
(and therefore) numerical accuracy decreases when the energy increases. 



Therefore a compromise should be achieved by selecting the levels in order to define the 
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objective function to minimize. 



We found that a reasonable compromise should be based on the following operations that 
differ with the selected rare gas: 

1. Pick some level and take the average of a number of higher ones. 

2. Pick two high levels for which the Quantum defect has stabilised within a given accu- 
racy. 

3. Pick a low level and a high one for which the Quantum defect has already stabilised. 

All the above operations should yield roughly the same values for the parameters before 
running the final check in order to test the accuracy of the obtained eigenvalues. 

The optimization program itself is based on a globally convergent method for solving 
non-linear system of equations: the multidimensional secant method developed by Broyden 
ril. It is based on a fast and accurate method for the iterative evaluation of the Jacobian 



of the objective function needed during the minimisation procedure. 

It is a Quasi-Newton method that consists of approximating the Jacobian and updating it 
with an iterative procedure. It converges superlinearly to the solution like all secant methods. 



There are several ways to perform the integration of the RSE on the basis of the 
CFM. One may use a fixed step scheme such as the explicit Runge-Kutta fourth order 
method (RK4) or a Variable Step with Control of Accuracy (VSCA). Optimization 
wise, the RK4 method is faster than VSCA but less accurate. To compromise, we 
first perform an initial search of the parameters with RK4 and then finalize the results 
using the VSCA method. Table | shows the parameters obtained for the first three Ryd- 
berg series of argon, optimized using the experimental energy levels tabulated by Moore [jl5| . 



To judge the accuracy of our optimization we compare in tables fT| and fT| our results 
for the energy levels and corresponding quantum defect obtained using the RK4 and VSCA 
integration scheme with those of Czydlik et al. |]T2[ . The results displayed in Tables || and 
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Series 


ei 


£2 


1=0 


3.625 


1.036 


1=1 


3.62 


1.06 


1=2 


3.6344 


1.036 



TABLE I: Czylik-Green parameters for the first three Rydberg series of Argon.. 



Experimental levels 


Calculated levels (RK4) 


Calculated levels (VSCA) 


-0.309522 


-0.214082 


-0.310563 


-0.124309 


-9.87422E-02 


-0.124506 


-6.76780E-02 


-5.69743E-02 


-6.76904E-02 


-4.25540E-02 


-3.70845E-02 


-4.25546E-02 


-2.92210E-02 


-2.60593E-02 


-2.92238E-02 


-2.13080E-02 


-1.93126E-02 


-2.13062E-02 


-1.62200E-02 


-1.48846E-02 


-1.62210E-02 


-1.27620E-02 


-1.18223E-02 


-1.27614E-02 


-1.03020E-02 


-9.61654E-03 


-1.03015E-02 


-8.49000E-03 


-7.97521E-03 


-8.49011E-03 


-7.11800E-03 


-6.72092E-03 


-7.11771E-03 



TABLE II: Comparison between the experimental and calculated energy levels of the Rydberg 
series of Argon. The levels calculated with the CFM are obtained either with fixed step (RK4) or 
variable step (VSCA) integration. All values in Rydbergs. 

PI clearly favor, as expected, the VSCA integration scheme. In contrast to the VSCA, the 
RK4 integration scheme is hmited to fourth order accuracy. Table |IIT| further shows the 
sensitivity of the QD to the numerical values of the calculated energy levels. Despite small 
differences between the energy values obtained with the RK4 and VSCA methods (Table 
the corresponding QD's largely differ. This is the main motivation for using QD rather 
than energy levels in our optimization procedure. 
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Experimental QD 


Calculated QD (RK4) 


Calculated QD (VSCA) 


0.202561 


0.838726 


0.205576 


0.163723 


0.817645 


0.165967 


0.156063 


0.810516 


0.156415 


0.152366 


0.807174 


0.152400 


0.150046 


0.805324 


0.150326 


0.149399 


0.804191 


0.149110 


0.148104 


0.803444 


0.148345 


0.148016 


0.802940 


0.147808 


0.147664 


0.802574 


0.147425 


0.147091 


0.802298 


0.147161 


0.147199 


0.802084 


0.146957 



TABLE III: Comparison between the experimental and calculated QD of the Rydberg series of 
Argon. The levels calculated with the CFM are obtained either with fixed step (RK4) or variable 
step (VSCA) integration. 

IV. SCATTERING METHOD AND RESULTS 

We now turn to the implementation of the optimized pseudo-potentials to the case of 
(e,2e) ionization in the equal energy sharing case (k^ = — k;,). Within a standard DWBA 
model and when a symmetric kinematic is considered (i.e. equal energy sharing between 
the two outgoing electrons) the exact unsymmetrized T-matrix element is approximated as 

{KK^Tmm}T^\^Ti:m^o) = (x-(ka)x-(k.)|V|</.L,MX+(ko)), (6) 

and both final state distorted waves, x~(ka) and x (kf,) are chosen as eigenfunctions of 
the electron-ion system. In this situation, the exact collision state with final state boundary 
condition, \l/*^~)(k(j, k^), is approximated as the product of two outgoing distorted waves. As 
for the calculation of the incoming distorted wave, the Furness-McCarthy local exchange 
approximation which, valid for scattering energies greater than 25eV 0], is commonly used 
is the calculation of the outgoing distorted waves. Furthermore, the distorting potential is 
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purely Coulombic asymptotically, and thereby neglects any e-e correlation in the final state. 

It is well known [0, that the product of two Coulomb waves satisfy the proper 
asymptotic condition provided the use of effective charges that satisfy the condition 

— h — — = \ I 1 . (7) 

ka kj, ka kfj jku kfti 

In the present case, we focus on the equal energy sharing geometry ka = — k;, and impose 
on the effective charge the condition 

z:" = z!" (8) 

Accordingly, we modify the form of our pseudo-potentials and replace in eq. (5), and for all 
r values, the unit charge by an effective charge ^'^■^■^ which satisfy (|^) and (^: 

Z'ff = -. (9) 
4 

This choice of screening charge insures that each outgoing electron experiences the exact 
classical force asymptotically. It was further suggested that the use of such effective charges 
suffices to take into account the electron correlation for the ionization of rare gas atoms in 



equal energy sharing geometry [IS 



In the spirit of a DWBA method based on pseudo potentials and to describe the ionization 
process of argon at energy beyond the validity of the Furness-McCarthy local exchange 
potential (below 25eV), we use the parametric form suggested by Jungen et al. ||2^ to 
describe the incoming distorted wave. This pseudo-potential is adjusted to reproduce the 
ab-initio elastic phase shift of McEachran et al. up to an energy of 30eV and includes 
static, exchange and polarization effects. At energy above 30eV, the Furness-McCarthy local 
exchange potential is used in the entrance channel. 

To gauge the improvement introduced using such a parametric potential to describe the 
incoming projectile, we show in Fig 2 a comparison of the elastic scattering phaseshift 
obtained using the pseudo-potential of Jungen et al, using the Furness-McCarthy local ex- 
change approximation and the ab-initio results of McEachran et al. which leads to scattering 
parameters in good agreement with the experimental values. The significant improvement 
introduced in the calculation of the elastic scattering phaseshift, and consequently in the 
determination of the radial part of the incoming distorted wave, strongly suggests that a 
superior description of the ionization process is to be expected. 
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FIG. 2: Comparison of the ab-initio elastic scattering phaseshift of McEachran et ai. (solid line) 
with the phaseshift obtained using the parametric potential of Jungen et al. (opaque symbols) and 
the Furness-McCarthy local exchange potential(filled symbols). 
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FIG. 3: Behavior of the TDCS as a function of scattering angle in equal energy sharing geometry 
for various values of the excess energy as indicated on each graph. DWBA calculations (thin line), 
DWBA calculations + pseudo-potentials (thick line), experiment (triangle). 
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Fig 3 shows the variation of triply differential cross section (TDCS) as a function of 
the scattering angle for various values of excess energy, equally shared between the two 
outgoing electrons. The results obtained using the pseudo-potentials defined in the previous 
sections for the incoming and outgoing distorted waves are compared to calculations where 
the Furness McCarthy local-exchange potential is used in both the entrance and the 
exit channels and to the experimental measurements available Since the experimental 
measurements are not absolute, they have been systematically normalized to the DWBA 
calculations at a scattering angle = 270 degrees except for an incident energy of Ei = 
17.8eV (i.e. Eexc = 2eV) and Ei = 19.8eV (i.e. Eexc = 4eV) where the experiment is 
normalized to obtain the best visual fit to the theory. 

For the highest incident energy considered here, Eq = 35.8 eV (i.e. E^^^^ = 20eV ), 
little differences in the magnitude of the TDCS is noticeable between the two theoretical 
models. In contrast, the use of the pseudo-potentials calculated in the previous section 
clearly improves the prediction of the shape of the TDCS. This is particularly the case 
of the angular region between 230 and 320 degrees. Although improved, the TDCS is 
however still slightly underestimated outside this particular angular range. We stress that 
at this particular energy, the Furness-McCarthy local exchange approximation is used in the 
calculation of the incoming distorted wave. Further, this result also validates our procedure 
for calculating the outgoing distorted wave for an energy interval (i.e 35.8eV for the incoming 
electron) where the Furness-McCarthy local exchange potential is still in its lower energy 
limit of validity. Similar remarks pertain for the magnitude of the TDCS at the next lower 
incident energy [E^xc = 8eV ). At this particular regime where the incoming distorted wave 
is now also described using pseudo-potentials, our optimized potentials largely improves the 
description of the experimental TDCS, especially in the angular region between 240 and 300 
degrees. 

At the two lowest energies considered here, marked differences are now noticeable be- 
tween the TDCS obtained using the two different theoretical models (i.e DWBA and 
DWBA-|-pseudo-potentials). First, the overall magnitude of the TDCS is largely reduced 
when parametric potentials are used. Additional calculations (not shown here) indicate that 
the use of effective charges is largely at the origin of this reduction. Unfortunately, the im- 
provement introduced can not be gauged against the experimental data as the measurements 
are not absolute. The reduction in amplitude of the variation of the TDCS where experimen- 
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tal measurements are available can be traced back to the influence of the pseudo-potential 



combined with the use of the effective charges in the exit channel ||2^. Finally, when the 
pseudo-potentials of Jungen et ai. are introduced to describe the entrance channel distorted 
wave little improvements are observed in the calculations of the TDCS where experimental 
measurements are available. Its effects on the calculated cross sections are mainly felt at 
the lowest and highest angular region where experimental measurements are not available. 



V. CONCLUSIONS 



In this work, we showed that QD information provides means to select and tune a class 
of pseudo-potentials. These can subsequently be used in ionization studies by electron 
impact with proper modifications to account for the electron correlation in the final state. 
Ionization of argon, in the equal energy sharing geometry, shows that this procedure provides 
an efficient treatment of the final state interactions in the energy region where the commonly 
used Furness-McCarthy local exchange potential is no longer valid, and yet allows for a 
much simpler alternative to the complete Hartree-Fock treatment. This might turn out to 
be particularly valuable when considering the ionization of molecular systems involving low 
energy electrons. We finally stress that the use of QD information, and consequently the 
application of the present method, would not be meaningful for energy above 30eV. 
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